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The synthesis of sol-gel silica materials doped with three different types of 
metallophthalocyanines has been studied. Homogeneous materials of good optical 
quality were prepared and the first optical limiting measurements of dyes in sol-gel hosts 
were carried out. The properties of these solid state limiters are similar to limiters based 
on phthalocyanine (Pc) in solution. 


Sol-gel silica materials containing copper, tin and germanium phthalocyanines 
were investigated. The initial step in all cases was to prepare silica sols by the sonogel 
method using tetramethoxy silane (TMOS), HC1 and distilled water. Thereafter, the 
synthesis depended upon the specific Pc and its solubility characteristics. Copper 
phthalocyanine tetrasulfonic acid tetra sodium salt (CuPc4S) is soluble in water and 
various doping levels (lx KHM to 1 x 10* 5 M) were added to the sol. The group IV Pc’s, 
(SnPc(OSi(n-hexyl)3)2 and GePc(OSi(n-hexyl)3)2, are insoluble in water and the process 
was changed accordingly. In these cases, the compounds were dissolved in THF and then 
added to the sol. The Pc concentration in the sol was 2 x 10' 5 M. The samples were then 
aged and dried in the standard method of making xerogel monoliths. 

Comparative nanosecond optical limiting experiments were performed on silica 
xerogels that were doped with the different metallophthalocyanines. The ratio of the net 
excited state absorption cross section (c e ) to the ground state cross section (a g ) is an 
important figure of merit that is used to characterize these materials. By this standard the 
SnPc sample exhibits the best limiting for the Pc doped sol-gel materials. Its cross 
section ratio of 19 compares favorably with the value of 22 that was measured in toluene. 

The GePc materials appear to not be as useful as those containing SnPc. The GePc doped 
solids exhibit a higher onset energy (2.5 mJ) and lower cross section ratio, 7. The 
CuPc4S sol-gel material has a still lower cross section ratio, 4, however, the 
tetrasulfonate groups make the dye soluble in water which greatly facilitates its 
incorporation into the sol-gel matrix. 
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The nonlinear transmission of CuPc4S in a pH 2 buffer solution and in a silica 
xerogel were compared. It is evident that the CuPc4S preserves its optical limiting 
behavior in the sol-gel matrix, indicating that the fundamental excited state absorption 
process is essentially the same for a molecule in solution or in the solid state. Although 
the spectroscopic details of energy level lifetimes are unknown, the significance is that 
passive optical limiting has been achieved in the solid state via incorporation of a dye into 
an inorganic host. The only compromise occurs at the extremely high energy regime 
where photobleaching is observed. This is a result of the limited mobility of the dye 
molecules in the solid silica host relative to a liquid host. The effects of 
photodegradation in the xerogel are additive, whereas the solution provides a supply of 
fresh molecules that are free to enter the active volume between pulses. 

In this program we supplied over 20 samples to JPL. The work was presented at 
the SPEE meeting in San Diego (July, 1992) and the paper was published in the 
conference proceedings: 

PD. Fuqua, K. Mansour, D. Alvarez, Jr. S.R. Marder, J.W. Peiry and 
B. Dunn, "Synthesis and Nonlinear Optical Properties of Sol-Gel Materials Containing 
Phthalocyanines," Proc. Sol-Gel Optics II. Vol. 1758, ed. by J.D. Mackenzie 
(SPIE, Bellingham, WA, 1992) pp. 499-506. 
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Introduction and Overview 


This report contains the results of research supported by the Jet Propulsion 
Laboratory (JPL) contract no. 959234, conducted by Prof. D. L. Mingori and 
K. F. Zimmermann. The report is in two parts. Part one is titled “System 
Identification Using Gram-Schmidt Orthonormalization”, and is concerned 
with the identification of linear, time-invariant systems, such as large flexible 
space structures. Part two is titled “Disturbance Rejection Using H 2 and 
H a o Control Design” , and is concerned with the design of time-invariant con- 
trollers for narrow and wide band disturbance rejection. During the period 
covered by this contract, the focus of the research was expanded beyond the 
original focus of system identification to include the subject of disturbance 
rejection. This change was initiated with the approval and support of the 
JPL Technical Manager, Dr. Mark Milman. 


Part I 

As stated above, part one of the report is titled “System Identification Using 
Gram-Schmidt Orthonormalization." System identification is the fitting of a 
chosen model structure to some data from the system. Some common models 
for linear, time-invariant systems are the Auto-Regressive Moving-Average 
(ARMA) model, the Auto- Regressive (AR) model and the Moving-Average 
(MA) model. This report focuses on these model types. 

In picking a model, the basis for the model must also be chosen. Often 
the basis that is chosen is not orthogonal. When the basis is not orthogonal, 
the model may be more sensitive to errors, and hence less “robust”. Part 
one of this report contains a general procedure to orthonormalize any basis 
when fitting a model to frequency domain data. 

The examples show that system identification using Gram-Schmidt or- 
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thonormalization has a great deal of potential. However, more work should 
be done in order to fully explore the usefulness of this procedure. 


Part II 

Part two of this report is titled “Disturbance Rejection Using Hi and 
Control Design.” To begin this section one needs to have a model of the 
system. This model could be found using the results of part one of this 
report. Then using Hi and Hoc control design techniques one can find a 
time-invariant controller that satisfies the design criteria. This report shows 
examples of both wide and narrow band disturbance rejection controllers. 

In this report we also explore the benefits of using H 2 and H^ control 
design techniques for multiple feedback control. It was found that for a given 
level of performance, single feedback controllers have higher gain than multi- 
ple feedback controllers. The lower gain of the multiple feedback controllers 
make them less sensitive to noise, and hence more “robust”. 


The authors would like to thank the Jet Propulsion Laboratory and in par- 
ticular Dr. Mark Milman and Dr. John Spanos for supporting this work. 
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Part I 
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Chapter 1 

System Identification Using 
Gram— Schmidt 
Orthonormalization 


Abstract 

Linear time-invariant plants are often modeled using Auto-Regressive Moving- 
Average (ARMA) models or Auto-Regressive (AR) models. Typically, these 
models are formed using basis functions corresponding to the z, s or 6 op- 
erator raised to powers. For example, a function A(p ) could be expressed in 
terms of a truncated power series: 


A(p) = a kP k 

k = 0 

where p represents either z, s or <5 operator, and the ajt’s are scalars. These 
bases may not be orthogonal. This is particularly apparent for s and 8. By 
replacing the basis {l,p, . . . ,p n } by an orthogonal basis, significant improve- 
ment can be achieved in the identification of a system. Presented in this 
report is a method of system identification using orthogonal bases found by 
Gram-Schmidt orthonormalization. 
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1.1 Introduction 


System Identification is the process of developing mathematical models using 
data from a system. The mathematical models are formed using combina- 
tions of functions. Two of the popular model types for linear time-invariant 
systems are the Auto-Regressive Moving-Average (ARMA) model and the 
Auto-Regressive (AR) model. A more complete discussion of these model 
types can be found in [8] chapter 4. This report focuses on using Gram- 
Schmidt orthonormalization in conjunction with these model types. 

Begin with the system, 

yip) = G(p)u(p) + H(p)e(p) (1.1) 

where y(jp), u{p ) and eip) are the output, input and the noise respectively; 
p represents either the z, s or 6 operator; G(p) and H(p) are the plant and 
noise transfer functions. Assume that y(p) and u{p) are known, while G(p), 
H(p) and e(p) are unknown. Consider a model in the form: 

yip) = G(jp , 9)u{p) + Hip , 9)eip, 0) (1.2) 

where G(p,0), H{p, 9) and eip, 9 ) are estimates of G(p), H{p) and e(p), 

and 0 is a vector containing the parameters of the model. Let eip, 9) = 
yip) - yip,0), and divide 1.2 through by H(p,0) to produce, 

y(p,e) = H(p,er'6(p,0)u(p) + (1 - H(p,er')y(p). (i.s) 

Define B(p,6 ) = Hip, 9)~ l G(p, 6), and Alp, 0) = Substitution 

using the definitions gives the following model, 

yip, 9) = B(jp, 9)u{p) + (1 - Aip,9)~ x )yip) (1.4) 

where the transfer functions A{p,9 ) and B(p,0) can be described by the 
following relationships, 


Aip,9) = 

1 + X >p* 

k=l 

n b 

(1.5) 

Bip,9) = 

£ bkP k - 
k=0 

(1.6) 


7 



For an ARMA model, n a > n& > 0, and for an AR model, n a > 0 and 
B(p,8) = 0. 

Notice that in the frequency domain p (i.e. z, s or S) is a function of the 
frequency u>, namely 

2 = e ju>T 
s = ju 
_ e jwT - 1 

6 = ~^f— 

From here on we will assume that p is a function of the frequency u (p = 
P(w))- 

Notice that the transfer functions A(p,&) and B(p,0) are represented by 
a series of functions. These functions, {p°,p l , . . . ,!?*}, form a basis for the 
series. Instead of using p to powers to form the basis of the series, a sequence 
of orthogonal functions could be used instead, such that 

A(p,0) = £a k j> k (uj) 

Jb=0 

B(p,0) = ^/3 k ^ k (u) 

fc= o 

where <t> k (u) is an orthogonal function. An orthogonal basis is considered to 
be a good choice for a basis, because it is typically less sensitive to errors, 
or more “robust” . In this report the terms “function” and “vector” are used 
somewhat interchangeably, since a discrete function (a function evaluated at 
several values of its argument) can be used to form a vector. 

This report discusses the properties of orthogonal functions, a method of 
generating orthogonal functions via Gram-Schmidt orthonormalization, and 
a method for using them in ARMA and AR models for the purpose of system 
identification. 


( 1 . 10 ) 

( 1 . 11 ) 


(1.7) 

( 1 . 8 ) 

(1.9) 


1.2 Theory 

1.2.1 Notation 

p Represents either the z, 6 or s operator. 
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z The Fourier series operator. 

6 The Delta transform operator. 

$ The Laplace transform operator. 
j The imaginary number, j = y/—l. 
uj Frequency in rad/ sec. 

T The sampling interval. 

0 A vector containing the parameters of the model. 
vj{x) A weighting function. 

w(x) A discrete weighting function (i.e. a vector containing the weights). 

W (x) A diagonal matrix containing the vector w(x) along the main diagonal. 

(j> n (x) The 71 th orthonormal function in an orthonormal function sequence. 

0„(x) The n th discrete orthonormal function (i.e. an orthonormal vector) in 
an orthogonal function sequence. 

T The transition matrix from the non-orthogonal basis to the orthonormal 
basis. 

1.2.2 Orthogonal Functions 

A function (or vector) sequence {0 O , <t >\ , 02, • • •} is said to be orthogonal if it 
obeys the following condition, 

\mtn (L12) 

where (•,•) is the inner product, and hm is a real finite non-zero scalar. If 
hm. = 1 then the sequence is called orthonormal. <pk is used to denote both 
4>k{x) and <£fc(x), where x is a real scalar and x is a real vector. 0*(x) is a 
discrete function produced by evaluating <pk(x) for x = Xi, where x, are the 
elements of x. 
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Continuous Functions 

For continuous functions, the inner product is defined as follows, 

(&»(*).&»(*)) = J a ™( x )<t>m( x )<f>n( x ) dx (1.13) 

where w(x) is a weighting function, and (•)* represents the complex conju- 
gate. The weighting function w(x) is real, non-negative and integrable on 
the interval (a, b), and w(x) > 0 on a sufficiently large subset of the interval 
such that the following is true, 

0 < [ w(x) dx < oo. 

J a 

Remark . In this treatment, only finite intervals will be dealt with (i.e. 
|a|,|6|<oo). 

Discrete Functions 

For discrete functions, the inner product is defined as follows, 

(4> m (x),0„(x)) = <f> m (x) H W(x)<j> n (x) (1.14) 

where 

X = [ Xi x 2 ■■■ x N ] r 

0m(5c) — ^ $m( x l) fim( x 2 ) 4 > m[ x N') j 

and W(x) is a diagonal matrix containing the weights, [-] w represents the 
complex conjugate transpose of a matrix. The weights, w(x) = diag(W(x)), 
are real, non-negative and 0 < ' 52^=1 w(xk) < x. 

1.2.3 Examples 

A couple of examples of common orthogonal functions are shown in order to 
familiarize the reader with these concepts. 
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The Exponential Function 

Let the orthogonal function be given by 4> n (w) = e jWr , the weighting func- 
tion by w(u) = 1, and the interval by (a, 6) = (0, ^). Then 



gjmu/T e — jnu/T fa 



(cos (muT) + j sin(muT))(cos(nuT) — j sin(nu/T)) du 
( cos(muT ) cos(nuT ) + sin (muT) sm(nuT )) du 


+j I (— cos (muT) sm(nujT) + sm(muT) cos (nuT)) du 
Jo 

2tr In 

/ cos((m — n)uT) du + j sin((m — n)ujT) du 
Jo Jo 


0 , m^n 

Y , m = n 


Hence the complex function e jnuT is orthogonal on the interval (0, y) using 
the weighting function w(u) = 1. 

The exponential function also forms a discrete orthogonal function when 
the frequencies are linearly spaced (i.e. u — ui = £ = 0, 1, . . . , N - 1), 

and the weighting is uniform (i.e. W(u>) = I, I is the identity matrix and 

u) = [ uiq u\ • ■ ■ us-\ ] )• This fact becomes important when perform- 
ing system identification in the z domain. A proof of this is shown in Ap- 
pendix 1.9.1. 


The Chebyshev Polynomials of the First Kind 

First the orthogonal function 0„(0) = cos (n9) with the weighting function 
w{9) = 1 on the interval (a, b) = (0,7r) will be developed as follows, 



cos (m9) cos (n9) d9 


1 f n 

- / (cos((m 4- n)9 ) 4- (cos((m — n)9) d9 

2 Jo 

0 , m n 

< 7 r , m = n = 0 . 

, f , rn = n ^ 0 


So, cos (n9) is an orthogonal function on the interval (0, tt) using a uniform 
weighting. Now, make the change of variable x = cos (9). The previous 
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development with this change of variables produces, 


J_{\-x 2 )-±T in {x)T n {x)dx=\ 


( 0 


IT 

n 

2 


, m 7^ n 
, m = n = 0 
, m = n 7^ 0 


where the weighting function is given by w(x) = (1 — x 2 )~ 5 and T n (x ) = 
cos(n cos -1 (x)) are the Chebyshev polynomials of the 1*‘ kind. In polynomial 
form, the Chebyshev polynomials are given by the following, 


7o (*) = 1 

Ti(x ) = x 

T 2 (x) = 2 x 2 — 1 

T$(x) = 4z 3 — 3x 

T n+ i(i) = 2 xT n (x)-T n . x {x). 


1.2.4 Gram-Schmidt Orthonormalization 

The Gram-Schmidt orthonormalization procedure is outlined as follows. Be- 
gin with a set of linearly independent functions (or vectors) {0 o, fa , , . . , 0 n }. 
Then form the corresponding orthonormal functions (or vectors) as follows: 



0 o 

00 

((00, 0o)) 5 

m 

0m+ 1 } (0m+l» 0fc)0fc 

Jfc=0 


0m+l = 


The resulting sequence {0o,0i, 


0m+l 

( (0m+l , 0m+l))^ 

. . . , 0 n } is orthonormal. 


1.2.5 The Weighted Least Squares Problem 

The system identification procedure in this report uses the weighted least 
squares solution. Some simple examples will be done in order to illustrate the 
characteristics of weighted legist squares solutions using orthogonal functions. 
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The Least Squares Problem 

Suppose it is desired to approximate a function /(x) using a finite orthogonal 
function sequence {(j>o(x) , fa (x) , . . . , <f> n (x)} . Let the following relationship 
exist, 

/(*) » Mfc(x). (1.15) 

k= 0 

By multiplying through by w(x)4>* m {x) (0 < m < n), and integrating, b m is 
found as, 

X? <£L(x)/(x)<ix 

b m = . (1.16) 

This corresponds to the weighted least squares solution to the problem, where 
the weighting is given by w(x). 

The Role of the Weighting Function w(x) 

The weighting function w(x) ends up being the weighting on the least squares 
problem of equation 1.15. In this section a couple of examples are done in 
order to illustrate the importance of the weighting function w(x). 

Example Using Chebyshev Polynomials of the 1 H Kind The poly- 
nomial f(x) = (x — l) 3 is approximated using the first three Chebyshev 
polynomials of the 1 st kind (i.e. T 0 (x) = 1, T\{x) = x and T 2 (x) = 2x 2 - 1) 
on the interval (a, b) = (—1, 1). The Chebyshev polynomials of the 1 st kind 
have a weighting function of w(x) = (1 — x 2 )~^ as shown in figure 1.1. The 
result is the following, 

/(*) ~ ~ z - ^(2x 2 - 1). 

Plots of the approximation and the approximation error are shown in figures 
1.2 and 1.3. 

Example Using Chebyshev Polynomials of the 2 nd Kind Again the 
polynomial /(x) = (x — l) 3 is approximated, but this time using the first 
three Chebyshev polynomials of the 2 nd kind (i.e. 7o(x) = 1, 7\(x) = 2x and 
T 2 (x) = 4x 2 — 1) on the interval (a, b) = (—1, 1). The Chebyshev polynomials 
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of the 2 nd kind have a weighting function of w(x ) = (1 — z 2 )2 as shown in 
figure 1.1. The result is the following, 


Weighting Functions of Chebyshev Polynomials of the 1st and 2nd Kind 



Figure 1.1: Plot of the Weighting Functions for Chebyshev Polynomials of 
the I st and 2 nd Kind 


f{x) « - 1 - + 7 -{ 2 x ) - ^(4x 2 - 1). 

Plots of the approximation and the approximation error are shown in figures 
1.2 and 1.3 for comparison with the results using Chebyshev polynomials of 
the 1** kind. 

Notice that the error in the middle part of the interval is less using the 
Chebyshev polynomials of the 2 nd kind (see figure 1.3). This can easily be 
explained by inspection of the corresponding weight functions (see figure 
1.1). Since the Chebyshev polynomials of the 2 nd kind place more relative 
weighting on the middle of the interval this method is more accurate there, 
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Function Value 


Approximating Functions Using Chebyshev Polynomials of the 1st and 2nd Kind 



Figure 1.2: Plot of the Approximating Functions Using Chebyshev Polyno- 
mials of the 1 st and 2 nd Kind 
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Errors in Approximations Using Chebyshev Polynomials of the 1st and 2nd Kind 



Figure 1.3: Plot of the Approximation Errors Using Chebyshev Polynomials 
of the l* 1 and 2 nd Kind 


16 



while the relative weighting on the middle of the interval is much less for the 
Chebyshev polynomials of the 1 st kind. 

1.3 Problem Statement 

The advantage of using orthogonal functions in system identification is that 
the approximation of the system data using a sequence of orthogonal func- 
tions is numerically more robust than an approximation using non-orthogonal 
functions. A paradigm of this is the approximation of a vector by a basis. 
Typically, the “best” basis to use is an orthogonal basis. Following is a 
development to illustrate where the difficulties arise. 

1.3.1 Fitting an Auto-Regressive Moving-Average Model 
to Frequency Domain Data 

An Auto-Regressive Moving-Average (ARMA) model in the s domain may 
be expressed as: 

n n 

ll(s, 0)=Yi t>kS k u(s) + ]T a k s k y(s ) 

k=0 k=\ 

where n a = rib = n for simplicity. Let s = ju:, and assume that u(ju) = 1. 

Then yijui) = G(uj), where Q(u) is the frequency domain data. If we wish to 
minimize the least squares cost function J = {y(jv) — y(jw,0)) H (y(ju)) — 
y(ju,& )) the following least squares problem must be solved, 

0 -[*‘ •*][ ~o\ 

where 0 is an appropriately sized zero matrix (a vector in this case), and 

= [0(«)]*6 

G{uo) 0 ••• 0 

0 Q{ u/,) ••• i 

; ••. o 

0 ••• 0 G(uh-i) 
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*6 

0 


' 1 jhjQ ■■■ (ju 0 ) n 

1 ju>i ■■■ (ju i) n 

. 1 juJN-l ••• (jUN-l) n . 

1 T 

ai a 2 • • • a n &o J . 


The regressor matrix contains the discrete basis functions for B(juj, 0). 
Let us take a closer look at these basis functions. It can be shown that 
column vectors which make up the regTessor matrix are not close to being 
orthogonal, and therefore, $b is not well conditioned. To illustrate this point 
a plot of the first few normalized column vectors is shown in figure 1.4. (In 



Figure 1.4: The First Few Column Vectors of the Matrix 

order to improve the numerics of the s domain formulation, the frequencies 
are often normalized such that ui is replaced by u> = If a function 


18 




which is orthogonal on the interval (u>o,u n-\) where used instead, the new 
regressor matrix 4»(, would be much better conditioned, and therefore the 
results would be numerically more robust. 

The first few Gram-Schmidt orthonormalized column vectors using the 
weighting w(u) = — 5 — r-, where u = u>i = — 1 4- 255 , i — 1, 2, • • • , 199, are 

shown in figure 1.5. This weighting corresponds to the weighting used by the 



Frequency 



Frequency 


Figure 1.5: The First Few Column Vectors of the Matrix 

Chebyshev polynomials of the l 5t kind. For the purpose of comparison, a plot 
of the magnitudes of the first few Gram-Schmidt orthonormal vectors and 
the first few normalized Chebyshev vectors are shown in figure 1.6. Notice 
the slight difference. The vectors formed using the Chebyshev polynomials 
of the 1 *‘ kind are not orthogonal, because we are taking u> at linearly spaced 
discrete points and the Chebyshev polynomials are not orthogonal as discrete 
functions for u linearly spaced. 
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Magnitude 


Comparison of Gram-Schmidt and Chebyshev Functions 

0.09 
0.08 
0.07 
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0.05 
0.04 
0.03 
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0.01 
0 

-1 - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 

Frequency 

Figure 1.6: Comparison of Gram-Schmidt (solid line) and Chebyshev Func- 
tions (dashed line) 
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1.4 


Background 

The use of particular orthogonal functions has been proposed by researchers 
in [6], [1]. Following is a brief discussion of other work which directly relates 
to the present work. 

1.4.1 The z Domain 

When performing system identification using frequency data in the z domain, 
the substitution z = e jujT is used. As shown in section 1.2.3, the functions 
1, e jwT ,. . . , e jnujT are orthogonal for the weighting w( u) — 1 on the interval 
(0,27r). Hence, the z domain already uses a sum of orthogonal functions 
to approximate the function. Although, it is only orthogonal when using 
uniform weighting and the complete range (0, 2ir) of linearly spaced points. 

1.4.2 The s Domain 

When performing system identification using frequency data in the s do- 
main, the substitution s = ju is used. In section 1.3.1, it was shown that 
the polynomials 1, ju ,. . . , (ju) n are not orthogonal. Dailey and Lukich [6] 
recommend the use of Chebyshev polynomials of the l j£ kind in this situa- 
tion. The following transfer function relationship between the polynomials 
of ju and the Chebyshev polynomials of u can be made, 


B(ju,0) 
A(ju , 0) 


£\,„ (-i) W + jiruo 

k even k odd 

£ n it=o {-l)*a k u k + jT, n k = 0 (-l) !i i 1 a k u k 

k even k odd 

H fc=o 0kTk(u)) + j 52 fc = o 0kTk{w) 
keven fc odd 

£ n fc=o a k T k (u) + j L n fc =o oc k T k {u) ' 

k even k odd 


Identification may be performed using the Chebyshev polynomials, and then 
converting the result from the Chebyshev domain to the s (or ju) domain. 

This method shows promise, but it may be possible to improve on it. 
Note that the Chebyshev polynomials of the 1 st kind are orthogonal on the 
continuous interval (—1, 1) with a weighting of w(u) = (1 — u; 2 ) - ^. Hence in 
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all cases of system identification the Chebyshev polynomials are only approx- 
imately orthogonal, because the data, Q(u >), is only known at a finite number 
a frequencies, and w(vk)T m (uk)T n (uk) i 1 0, for all m ± n. Also, if a 
weighting other than w( u) = (1 — is used the Chebyshev polynomials 

of the l 4t kind are not orthogonal. 

1.4.3 The 6 Domain 

Recall that 6 = £ jr-. Therefore in the frequency domain, 6 = ? - y 1 . Bayard 
[1], recommends the use of the Chebyshev polynomials of the 1 st kind in the 
6 domain when the following condition is met: The sampling rate is very fast 
(or T very small) when compared to the frequency of the highest mode of 
the system. When this condition is met the 6 domain is very similar to the s 
domain, hence the Chebyshev polynomials are approximately orthogonal in 
the 6 domain. 

The Gram-Schmidt procedure outlined in this report allows one to find 
an orthonormalization of any domain, and with any weighting. 


1.5 System Identification Using Gram— Schmidt 
O rt honor malizat ion 

The first step in system identification is to obtain data, and in this procedure 
we begin with frequency data. The procedure of obtaining data falls under 
the heading “experiment design”. Ljung [8] devotes a chapter to this subject. 

Here we will briefly mention some of the more popular methods of obtaining 
frequency domain data. 

Perhaps the most obvious method is a sine dwell or sine sweep method. 

For the sine dwell method, the system is forced with a sine wave at constant 
frequency and magnitude until the transients die out, and then the phase 
and magnitude of the output are found. This is then done for additional 
frequencies until an adequate frequency record is found. The sine sweep 
method is similar except that the forcing frequency is not constant, but 
adjusted slowly. 

Another method is to use time domain data from the input and output. 

This data is transformed into the frequency domain using the Discrete Fourier 
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TYansform (DFT), and the output frequency data is divided by the input 
frequency data in order to find the system frequency data. Care must be 
taken to insure that the input used is “rich” enough in the frequency domain 
so that the results are accurate. Also, There are several other issues that 
must be dealt with here. Two of the most important issues are “aliasing” 
which is related to the high frequency modes of the system and sampling 
rate, and “windowing” which is important if the input data record does not 
contain an integer number of periods of the input signal. 

1.5.1 The General Procedure 

1. Obtain frequency data G(u). 

2. Choose the Model Type. We will assume the most general model here, 
the ARM A model, 

yip, 9) = - 2; a*p‘SM + £ hr/-. (1.17) 

1 fc =0 

where we assume that u(p) = 1 and y(p) = G(co)- In order to keep 
this development general, the frequency dependence of p will be shown 
using the notation p = p( u). For simplicity, let n a = rib = n. This leads 
to the following equation for the weighted linear least squares problem, 

0 = W(u;)*[* o * 6 ] * (1.18) 

where 0 is an appropriately sized matrix of zeros (in this case a column 
vector of size (N — 1) x 1), 

id(cjo) 2 0 ••• 0 

0 w(uj\)* : 

: ••• ••• 0 

0 ••• 0 iu(u;jv-i) 5 

GM G{vo)p{uo) ■■■ G{vo)p{uJo) n 

G(uj i) G(uji)pM ••• G(ui)p{ui) n 

= . . 

. G(u>N-l) G(u>N-\)p(^>N-i) ••• G(<^N-l)p(^N-l) n 
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' 1 p(u o) • • • p(u o) n 

1 p(u i) ••• p(cui) n 

?6 = . 

.1 p(u N - 1) ••• p(u N -i) n m 

9 — ^ ®1 ^2 ' ' * On ^0 ^1 ‘ ‘ ‘ j 


3. Orthonormalize the Regressor Matrix. At this point we have a few 
choices of how to proceed. Orthogonalization could be done with re- 
spect to W(u)) 2 $ 0 , W(cj)^$ 6 or W(o>)5 [ $ a ]• We will choose 
the following method: The orthonomalization will be done with respect 
to W(u>)i& b such that $ b = $ b T and 9(*W(uj)9b = I. Therefore, 
our least squares problem can be written as follows, 


0 


[ 9 a T * b T ] 


1 

T -1 

0 


' -1 ■ 

i 

0 

T -1 


9 


W{ w)* [9a Qb 


c 

0 


(1.19) 

(1.20) 


where $ b is orthonormal, $ a is not orthonormal, c is a known con- 
stant (because T -1 is a lower triangular matrix), and 6 contains the 
parameters of the orthonormalized model. 


4. Solve the Linear Least Squares Problem. Using Matlab notation, the 
weighted least squares equation is the following, 

cW(«)4,(:,l) = lV(w)l [ * 0 (:,2:n 0 ) * b ] 9 (1.21) 

where $„(:, 1) is the first column of $ a and $ o 0,2 : n a ) is the rest of 

* a . 

5. Convert the Model Back into the Original Domain (find 9). The model 
parameters in the original domain can easily be found as, 


9 = 


T(2:n+l,2:n+l) 0 

0 T 


9. 


( 1 . 22 ) 
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Remark 1: Although 9 is real, 9 may not be real. In order to insure that 
T is real, and hence 0 is real, the frequency data Q{uj) and the weighting 
w(u)) should be chosen to be symmetric about u = 0. 

Remark 2: It should be pointed out that orthonormalizing with respect 
to essentially finds an orthonormal basis made up of polynomials 

of p(u), i.e. 

t 

M“) = a kp k (u ) 

k = 0 

where a* is a scalar (and real if Remark 1 is followed). 


1.6 Examples 

The following examples use a single-input single-output state space model 
generated from a finite element model of the Jet Propulsion Laboratory 
(JPL), Control and Structure Interaction (CSI), phase B testbed. The model 
was obtained courtesy of JPL’s CSI group, and in particular Dr. Sam Sirlin. 
The model has 33 detectable modes, in the frequency range from 0 to 128Hz. 

For all the tests, the assumed order of the model began at 33 and was 
decreased until a stable minimum order model was identified. By “minimum 
order” we mean a model which does not contain any pole zero cancellations. 
Also, for all the tests 513 linearly spaced frequencies points in the range 
(~ f>f) symmetrically positioned about u = 0 were used (i.e. u = ut = 
(-7r + %j£)/T, t = 0, 1, . . . ,512). Once A(p, 9) and B{p,9) are found, they 
are used to form the model 

^ 9) =Wi u{r) - (l - 23) 

This model is then compared to the original data. 

1.6.1 z Domain Identification 

The first identification was done using a uniform weighting (WCu;) = I). 
Both the z domain and orthonormalized identified models are 33 rd order. 
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The results are shown in figure 1.7. The magnitude of the system data 
||^(cj )||2 is shown as well as the errors in the two identified models. The 
errors are defined as G(u)) — y(e :LJT , 0). Notice that the identified models are 


Frequency Response of ID in the z and Orthogonal Domains (W=I)) 



Figure 1.7: Frequency response magnitude plot of the plant, the error in the 
z domain identified model and the error in the orthonormalized identified 
model. The identification used uniform weighting. 

almost indistinguishable. This result should be expected since the z domain 
is orthogonal for unit weighting (although it is not exactly orthogonal here 
since we are not using the proper u ; in order to make it orthogonal). 

Next, we will use the weighting w(u;) = ||C?(u>)|| 2 . This is a reasonable 
weighting to use, because the weighting of w(u) = |j ^ g \\l gives us an 

unbiased model, but since we do not know A(e jU)T , 0) we use £(u>) in its 
place. The models are both 32 nd order. The results are shown in figure 1.8. 
The difference between the models is noticeable, but it is not apparent that 
one model is superior to the other. 
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Frequency Response of ID in the z and Orthogonal Domains (W=IGI)) 



Figure 1.8: Frequency response magnitude plot of the plant, the error in the 
z domain identified model and the error in the orthonormalized identified 
model. The identification used a weighting of ||£(^)|| 2 - 
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1.6.2 s Domain Identification 


All the identifications in this section use unit weighting, and the frequencies 
are the same as those in the z domain identification. 

First, identification was done using the s domain and the Gram-Schmidt 
orthonormalized s domain. The s domain model is 3 rd order, and the or- 
thonormalized s domain model is 28 th order. As can be seen in figure 1.9, 
the s domain identification model is not good at all, while the orthonormal- 
ized s domain model is much better. 


Frequency Response of ID in the s and Orthogonal Domains (W=I) 



Figure 1.9: Frequency response magnitude plot of the plant, the s domain 
identified model and the orthonormalized s domain identified model. The 
identification used uniform weighting. 

Next, identification was done using the scaled s domain and the Gram- 
Schmidt orthonormalized seeded s domain. The scaled s domain is used 
to improve the numerics of the s domain formulation. Scaling is done by 
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replacing u; by u = — — , where Umax = max(\u\). The scaled s domain 
model is 24 th order, and the orthonormalized scaled s domain model is 28 th 
order. As can be seen in figure 1.10, the scaled s domain identification model 
is much improved over the unsealed s domain, but still not as good as the 
orthonormalized scaled s domain model, or the orthonormalized s domain 
model. 


Frequency Response of ID in the Scaled s and Orthogonal Domains (W=I) 



Figure 1.10: Frequency response magnitude plot of the plant, the scaled s 
domain identified model and the orthonormalized scaled s domain identified 
model. The identification used uniform weighting. 

Figure 1.11 displays the error in the orthogonalized unsealed and scaled 
s domains for comparison. Both models are 28 th order, and very similar. 
This result is not unexpected since the Gram-Schmidt orthonormalization 
procedure does its own scaling. Prescaling seems to have little effect on this 
method. 
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Frequency Response of ID in the s Orthogonal and Scaled s Orthogonal Domains (W=I) 



Figure 1.11: Frequency response magnitude plot of the plant, and the errors 
in the orthonormalized s domain identified model and the orthonormalized 
scaled s domain identified model. The identification used uniform weighting. 
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1.6.3 5 Domain Identification 

All the identifications in this section use unit weighting, and the frequencies 
are the same as those in the z domain identification. 

First, identification was done using the <5 domain and the Gram-Schmidt 
orthonormalized 6 domain. The highest order stable 6 domain model is 4 th 
order. As can be expected from a 4 th order model approximating a 33 rd 
order plant, the result was very poor. The orthonormalized 6 domain did 
not produce a stable model for any order from 1 to 33. For the higher order 
attempts, the orthonormalized 6 domain models fit the data well, but were 
unstable. 

Next, identification was done using the scaled 6 domain and the Gram- 
Schmidt orthonormalized scaled 6 domain. Like the scaled s domain, the 
scaled 6 domain is used to improve the numerics of the domain. Scaling 
is done by replacing v = e, ^~ l by 0 = , where ^ = max(\v\). The 

scaled 6 domain model is 23 rd order, and the orthonormalized scaled 6 domain 
model is 22 nd order. The results are shown in figure 1.12. The scaled 6 domain 
identification model is much improved over the unsealed S domain, and also 
better than the orthonormalized scaled S domain model. 

The Gram-Schmidt orthonormalization does not seem to be of benefit for 
this type of model in the 6 domain. This result is surprising considering the 
success in the s domain. Part of the problem is that the higher order, better 
fitting models in the orthogonalized domain were unstable. So now we will 
compare the method once again, but this time allow unstable models. All 
the models are 33 rd order and unstable. 

Figure 1.13 displays the 6 domain and the Gram-Schmidt orthonormal- 
ized 6 domain models. The Gram-Schmidt orthonormalized S domain model 
fits the data much better than the <5 domain model. Figure 1.14 displays 
the scaled 6 domain model and the Gram-Schmidt orthonormalized scaled 6 
domain model. These models are very similar. For a better comparison, the 
model errors are shown in figure 1.15. The orthonormalized scaled 6 domain 
model is better in the low frequencies, but worse in the high frequencies com- 
pared to the scaled 6 domain model. When the model order was permitted 
to increase the errors were much less, but the models were unstable. 

Although system identification using Gram-Schmidt orthonormalization 
did not perform well for this model in the 6 domain, it should be studied 
further. We will continue to research this problem with the intent of finding 
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Frequency Response of ID in the Scaled delta and Orthogonal Domains (W=I) 



Figure 1.12: Frequency response magnitude plot of the plant, the scaled 6 
domain identified model and the orthonormalized scaled 6 domain identified 
model. The identification used uniform weighting. 
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Frequency Response of ID in the delta and Orthogonal Domains (W=I) 



Figure 1.13: Frequency response magnitude plot of the plant, the 6 domain 
identified model and the orthonormalized <5 domain identified model, both 
33 rd order and unstable. The identification used uniform weighting. 
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Frequency Response of ID in the Scaled delta and Orthogonal Domains (W=I) 



Figure 1.14: Frequency response magnitude plot of the plant, the scaled <5 
domain identified model and the orthonormalized scaled 6 domain identi- 
fied model, both 33 rd order and unstable. The identification used uniform 
weighting. 
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Frequency Response of ID in Che Scaled delta and Orthogonal Domains (W=I) 



Figure 1.15: Frequency response magnitude plot of the plant, and the errors 
in the scaled 6 domain identified model and the orthonormalized scaled 6 
domain identified model, both 33 rd order and unstable. The identification 
used uniform weighting. 
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the reasons for the difficulties and solutions to them. 


1.7 Multiple-Input Multiple-Output Gram— 
Schmidt Orthonormal Identification 

The MIMO system is the following, 

y (p) = G(p)u(p) + H (p)e(p) (1.24) 

where y (p), u(p) and e(p ) are column vectors containing the output, input 
and noise respectively; G (p) and H(p) are the plant and noise matrix trans- 
fer functions. As in the single-input single-output case, we begin with the 
following model, 

y(p) = G(p,0)u(p) + H(p, 0)e(p, 0) (1.25) 

where G(p, 0), H (p, 0), and e(p, 0) are estimates of G(p), H(p), and e(p). 
Let G(p,6) - and assume that H(p, 0) = , where A(p,0) is 

a scalar transfer function, B (p,0) is a matrix transfer function, and I is an 
appropriately sized identity matrix. The assumption that H(p, 0) — A{ ^q^ 

greatly simplifies the result. Now, by multiplying through by A(p, 0), and 
substituting e(p,0) = y(p) — y(p, 0) the following equation is found 

y(p,0) = B(p,0)u(p) + (1 - A(p,0)~ l )y(p) (1.26) 

where 

B n (p,0) B x2 (p,0) ••• B lni (p,0) 

B (p,0)= t (127) 

. B nol (p,0) ••• B noni (p,0) 

rii is the number of inputs and n a is the number of outputs. A(p,0) and 
Boi{p,0) can be described by the following relationships, 


A(p,0) 

= i 

*:=i 

Kb 

(1.28) 

Bo%(p,0) 

= Y,boi,kV k 

k= 0 

(1.29) 
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where i = 1 , 2, . . . , n» and o = 1 , 2, . . . , n 0 . To simplify the notation, we have 
made all the series which approximate the transfer functions the 

same length, n 

In order to set up the least squares problem, we need n* x n a frequency 
data records £»(uj), where Goiiu) is the frequency data from input i to output 
o. Setting up the weighted least squares problem associated with 1.18 results 
in the following equation, 


0 = W(u>)$$ 


(1.30) 


where W(cu) is size ( nin 0 N x Tiin 0 N), and is size (nin 0 N x ( n a + ribn 0 )) 



‘ W( u)h 

0 

0 ■ 


W ( U >)5 = 

0 

W( u>)5 ••• : 

0 



0 


0 W(u)l J 


$ = 

r [*.]» 
[*.] a 

i 



. (*«L 

l*>). 


[*.]. = 

a, o2 




. &a,orii . 



. = 
^ 0,01 — 

GoiM 

Goi(w l) 

GoxMpM 

GoxMp{u>i) 

Goi(vo)pM na 

Gox(ui)pM na 


|_ Goii^N- 

-l) Goi(uN-l)p(UN-l) ' 

•• Goi(v N - l)p(u> N -i) na . 


r** o 


0 1 


0 ••• 0 $ b 
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1 p(wo) • • • p(u/ 0 ) nt 

1 p(wi) ••• p(uq) nt 


L 1 p(u N -i) ■■■ p(u N -i) nb J 


a 

bn 

bn 



a 



In the MIMO problem just formulated, it seems to make the most sense to 
orthonormalize with respect to This would produce the following 

equations, 


0 = W(w)5 


W, 1*6] 

[*»] 


* 


°J2 


[T] [T 


-n 


0 = W(w)J 


[*a] n „ [*>] 

*a| 


$(>! 

4>(,i 




c 

0 


l*«L (*•) 


0 = W (u)l*0 


c 

9 


-1 

0 
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where 


••• 01 


0 

0 T 

[$ 0 j is not orthogonal and [$*,] is orthonormal. So, like the SISO case, we 
need to solve the least squares problem 

c*.(:, 1) =[*.(:, 2: n.) ] 0 , (1.31) 

for 0 and then convert to the original domain using [T]. 


[T] = 


T 0 
0 T 

0 ••• 




1.8 Conclusions for Part I 

This report outlines a general method for doing system identification of linear 
time-invariant systems based on Gram-Schmidt orthonormalization. The 
advantage of this method is that it allows one to orthonormalize in any 
domain with any weighting. 

System identification using Gram-Schmidt orthonormalization achieved 
significant improvement in the identification of models in the s domain, as 
compared to the more standard methods. The use of Gram-Schmidt or- 
thonormalization in the z domain showed no improvement, and in the 6 
domain we encountered problems in applying the method. More research is 
necessary to fully evaluate the benefits and difficulties associated with this 
method of system identification. 


1.9 Appendix 

1.9.1 Orthogonality of Discrete Exponential Function 

Claim: The discrete function e jU}T , where u> = [ uo uq •• • u>n-i ] and 
ue = is orthogonal with unit weighting, w(ue) =1. 
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Proof: The discrete function e jU,T is orthogonal with unit 
the following holds, 

weighting if 

C _ ( JUTmsH 0 jU)Tk _ J , k = TO 

- (e ^ } 6 “ { 0 ,k?m- 

(1.32) 

where hm is a real scalar and 0 < hm < oo. Changing from vector notation 

to summation notation, 


N- 1 


Sn* = Je- iU(T V“< Tt 

(1.33) 

<=0 


N-l 


= £ e iu ' <7 ’ 9 

*=0 

(1.34) 


where q = k — m. At this point, it is obvious that if k = m (q = 0), then 
Smk — N+l ( hm = N+ 1). Now, we must show that if k ^ m, then S mk = 0. 
This can be shown as follows: Expand 1 .34 

= 1 + e j “ lTq 4 • • • 4 4- e juN ~ lTq , 

then multiply by e juiTq , 

e^ xTq Smk = + gi("l+wi )Tq _j_ . . . _|_ gX^i+w/v -i)Tq gj(wi+uw -i)Tg 

Notice that + ue = u>e+u and that = u 0 = 0 rad, so that 

e Mr *S*»* = e juiTq 4 e junTq + ••• 4 e i “ N - lTq 4 1 = S mk - 

Since e jwiTq Smk = S m k , and e ju>lTq ^ 1 for q ± 0, this implies that S m k = 0 
for q ^ 0. 
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Chapter 2 

Disturbance Rejection Using 
H .2 and Hoc Control Design 


Abstract 

Wide and narrow band disturbance rejection is explored using H 2 and Hoc 
control design. This design method allows for multiple feedback design, there- 
fore, the advantages of using primary and secondary feedback measurements 
verses only primary feedback measurements for disturbance rejection are also 
discussed. The primary measurement is a measurement of the signal which 
is to be controlled, and the secondary measurement is a measurement which 
is closely correlated to the disturbance. 


2.1 Introduction 

This report discusses research disturbance rejection using H 2 and Hoo control 
design. First, the design procedure is developed. In this procedure both the 
inputs and the outputs of the closed loop system are weighted. Weighting the 
inputs corresponds to including a priori information in the system for control 
design. Weighting the outputs corresponds to enforcing the performance 
criteria. Examples are then used to illustrated the method. 

Also, the use of primary and secondary feedback measurements verses 
only primary feedback measurements is explored. Primary measurements 
correspond to outputs one wishes to control, while secondary measurements 
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correspond to outputs one does not care to control but which are highly 
correlated to the disturbance. The motivation for using secondary measure- 
ments is that the more information the controller has, the better the control 
should be. For instance, in disturbance rejection one obvious measurement 
would be the output of the sensor corresponding to the quantity which is be- 
ing minimized. This is called the primary output. This may be a force sensor 
at the location where it is desired to have zero force. Another output which 
is desired is the disturbance itself. This is not usually a practical measure- 
ment, but perhaps it is possible to place a sensor “near by” the disturbance, 
so that this output is well correlated with the disturbance. This is called the 
secondary output. An example is included which illustrates the advantage 
of using a primary (controlled) output and a secondary (disturbance) output 
as feedback. 


2.2 General Discussion 

First, notation and the problem framework are defined. 


2.2.1 Notation and Problem Framework 


Begin by describing the system in figure 2. 1 , where P is the open loop system, 
K is the controller, w is the exogenous input, u is the controlled input, z 
is the controlled output, and y is the sensor output which is fed back to 
the controller. P and K are matrix transfer functions, while w, u, z and 
y are vectors. Note that in general w contains all exogenous inputs (i.e. 
reference inputs, disturbances, sensor noise, process noise), and z contains 
all controlled outputs for the performance criteria. Therefore, z may include 
the controlled input tt if one of the design criteria is to limit u in some way. 
P can be written as follows, 


P = 


Pu 
. P2 ‘ 


Pn 
P 22 


where P n corresponds to the open loop transfer function matrix from w to 
2 , and so forth. Thus, the system equations are the following, 


Z = P\\W + P l2 u 
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w 


z 


u 



y 


Figure 2.1: Closed Loop System 

y = P 21 W + P 22 U 
u = Ky. 

The closed loop transfer function matrix T from the exogenous input w to 
the controlled output z is given by, 

T=P n + P n K(I - P22K)- l P 2l . (2.1) 

The H 2 control design finds the stabilizing controller K which minimizes 
the 2-norm of T (i.e. minK\\T\\ 2 ). 

The Hoc control design finds a stabilizing controller K such that the oo- 
norm of T is less than or equal to one (i.e. j|T’|| 00 < 1). 

In order to include a priori information concerning the exogenous inputs 
and to enforce the design criteria, input and output weightings are added to 
the system. 

2.2.2 Input and Output Weightings 

Weightings can be added to the system shown in figure 2.1. These weightings 
axe transfer functions, and can be used to incorporate a priori knowledge into 
the system and to enforce design criteria. The placement of the weightings 
is shown in figure 2.2, where Wi is the input weighting matrix and W a is the 
output weighting matrix. Wi and W a are diagonal transfer function matrices. 
Now, the system equations are the following, 

W~ l z = P n W i w + P l2 u 
y = P 2 iWiW + P 22 U 
u = Ky. 
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Figure 2.2: Closed Loop System with Weightings 

Therefore, the weighted closed loop transfer function matrix T from the 
disturbance w to the weighted controlled output z is given by, 

T=W 0 (P n + P 12 K(1 - P 22 K)~ l P 2X ) Wi = W 0 TWi. (2.2) 

Following is a discussion of the roles of the input and output weightings. 

Input Weighting Matrix 

The input weighting matrix is used to incorporate a priori knowledge of the 
inputs into the system. For instance, if we know that the sensor noise is below 
some value, then we would use that value as our weighting for the input which 
corresponds to the sensor noise. For example, if sensor noise is given as in 
figure 2.3, then a transfer function is used which closely resembles this as the 
weighting on the sensor noise input, in order to include this information in 
the design. The reason to use these weightings is that the H 2 and H designs 
work to produce a controller K such that T is minimized in some respect. So, 
in a rough sense the H 2 and Hoc designs assume that the exogenous inputs 
w are full spectrum with a magnitude of one, so that the outputs are given 
by the transfer functions in T. Therefore, the input weightings are used to 
incorporate what we know about the inputs into the system. Examples are 
given in section 2.3 of this report. 

Output Weighting Matrix 

The output weighting matrix would be used to enforce design criteria on the 
closed loop system. For instance, in the case of disturbance rejection where 
the system output z is the plant output, and it is desired to minimize the 
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Figure 2.3: Example of Sensor Noise Spectrum 

output at high frequencies, but the system can tolerate larger outputs at low 
frequencies, then the weighting shown in figure 2.4 may be appropriate. This 



Figure 2.4: Example of Output Weighting Spectrum 

weighting forces the unweighted closed loop transfer function T to be small 
where W a is large, because the control design is done on the weighted system 
T = W a TWi. Therefore, weighting the output corresponds to enforcing the 
design criteria. 
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2.3 Examples 

The test model is a one mass model with one control input and two sensors. 
The model is developed and then used to illustrate the design method. The 
designs are for rejection of two types of noise, wide band and narrow band 
disturbances. 

2.3.1 Test Model 

The one mass model is shown in figure 2.5. The constants used in the model 



Figure 2.5: One Mass Model 

are given in table 2.1. The damping coefficients correspond to 2% damping 
of the modes, and the outputs (yi, 1 / 2 ) correspond to the force at the wall 
and the acceleration of the mass. The sampling interval is T = 0.001. The 

m 1 = 2.18 kg 
ki = 14 N/mm 
ci = 6.9880 kg/s 

Table 2.1: Values Used in One Mass Model 

block diagram for this system is shown in figure 2.6. The dependence of the 
transfer function on s, the Laplace variable, is not stated explicitly to keep 
the notation cleaner. A bode plot of the plant transfer function from the 
disturbance d to the controlled output yi is shown in figure 2.7. 
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The complete closed loop system block diagram is shown in figure 2.8. 
This system includes the input and output weightings, a delay in the response 



Figure 2.8: Closed Loop System Including Weightings 


of the controller, and sensor noise. It does not include process noise at this 
point. The controlled outputs are y\ and u. So, corresponding to the earlier 
system description P, K,w, w, u, i, 2 and y are given as follows 
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where D(T d ) is a delay of time T d , and W x and W Q are diagonal transfer 
function matrices containing the appropriate weightings. By absorbing the 
weightings and the delays into the open loop plant P, a new open loop system 
P is formed given by the following 
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2.3.2 #2 and #<» Control Design Tests 

Two different tests axe run. The first is a demonstration of wide band dis- 
turbance rejection, and the second is a demonstration of narrow band dis- 
turbance rejection. The tests are run using the one mass model described 
previously. And, the delay, D(Td ), is approximated using a second order 
Pade approximation, where Td = 0.001. The tests are run using the Matlab 
Robust-Control Toolbox [3]. 

Wide Band Disturbance Rejection Using H 2 and //,» Control design 

H 2 Control Design for Wide Band Disturbance Rejection The first 
step is to set up the input and output weightings. These are shown in figure 
2.9. The input weightings were all chosen to be constant transfer functions 
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H2 Design: Wide Band Disturbance Rejection: Input Weightings 
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Figure 2.9: Input and Output Weightings for H 2 Design Wide Band Distur- 
bance Rejection 
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(no dependence on frequency). The weighting on the disturbance cf is 1 , 
while the weighting on the two sensor noises n sj and n a2 is equal and is 0.001. 
This corresponds to most of the noise in the system being attributed to the 
disturbance, and not the sensor noise. The output weightings correspond 
to a design criteria which demands better disturbance rejection at higher 
frequencies, and less control force at higher frequencies. 

The closed loop and open loop frequency response is shown for comparison 
in figure 2.10. At high frequencies the disturbance rejection is quite good, 


H2 Design: Wide Band Disturbance Rejection 



Figure 2.10: Comparison of Open Loop and Closed Loop Transfer Functions 
for H 2 Design Wide Band Disturbance Rejection 

and not as good in low frequencies. The poor results at low frequencies can 
be explained by examining the model (figure 2.5). For a constant disturbance 
force d, the control force u would have to be infinite in order to have yi equal 
to zero. Since the control force u is part of the performance measure, by 
being included in z, solutions which rely on unrealistic control forces Eire 
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ruled out. The closed loop transfer functions from the sensor noises to the 
output j/i were inspected to insure that they contribute very little to the 
output yi when compared to the effect of the disturbance d. 

A frequency response plot of the controller is shown in figure 2.11. 



Frequency (Hz) 


Figure 2.11: Controller Transfer Function for H 2 Design Wide Band Distur- 
bance Rejection 


Hoc Control Design for Wide Band Disturbance Rejection Again, 
the first step is to set up the input and output weightings. The same input 
weightings are used from the H 2 design case, but the output weightings are 
scaled so that a design solution can be found. The seeded output weights 
are shown in figure 2.12. It was possible to find.a solution without changing 
the weighting on the controlled input u, but the weighting on y x had to 
be so relaxed that the solution was to use a very low gain controller. This 
corresponds to essentially disconnecting the controller. Also, the weighting 
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H-inf Design: Wide Band Disturbance Rejection: Output Weighting 



Figure 2.12: Output Weightings for Design Wide Band Disturbance 
Rejection 
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on y i could be raised, but at the expense of having to decrease the controlled 
input weighting, which corresponds to allowing an even greater control force. 
The weighting used for W o2 , the weighting on u, allows the design to have 
a much greater controlled input than the H 2 design. This is not desirable. 

The closed loop and open loop frequency response is shown for comparison 
in figure 2.13. 


H-inf Design: Wide Band Disturbance Rejection 



Figure 2.13: Comparison of Open Loop and Closed Loop Transfer Functions 
for Hoc Design Wide Band Disturbance Rejection 

A frequency response plot of the controller is shown in figure 2.14. 

Narrow Band Disturbance Rejection Using H 2 and H «, Control 
design 

H 2 Control Design for Narrow Band Disturbance Rejection The 
input and output weightings are the same as those for the H 2 design for wide 
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Figure 2.14: Controller Transfer Function for H x Design Wide Band Distur- 
bance Rejection 
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band disturbance rejection except for the input weighting on the disturbance 
d. The input weighting on the disturbance d is formed to account for a 
narrow band disturbance at 10 Hz, and broad band low level noise. The 
weightings are shown in figure 2.15. 
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H2 Design: Narrow Band Disturbance Rejection: Input Weightings 
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Figure 2.15: Input and Output Weightings for H 2 Design Narrow Band Dis- 
turbance Rejection 

The closed loop and open loop frequency responses are shown for com- 
parison in figure 2.16. 

A frequency response plot of the controller is shown in figure 2.17. The 
weighting of the disturbance causes the controller to have an inverted notch at 
10 Hz. This is difficult to see on the magnitude plot, but is more apparent on 
the phase plot. This result can be explained by the Internal Model Principle 
[9], where the dynamics of the disturbance is included in the controller, so 
that in the closed loop system the dynamics of the disturbance are canceled 
by those of the controller. 
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Figure 2.16: Comparison of Open Loop and Closed Loop Transfer Functions 
for // 2 Design Narrow Band Disturbance Rejection 
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Figure 2.17: Controller Transfer Function for H? Design Narrow Band Dis- 
turbance Rejection 
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Hoo Control Design for Narrow Band Disturbance Rejection The 
input weightings are the same as those for the H 2 design for narrow band 
disturbance rejection. Again, the output weightings are scaled so that a 
design solution can be found. The scaled output weights are shown in figure 
2.18. As with the wide band case, it was possible to find a solution without 


H-inf Design: Narrow Band Disturbance Rejection: Output Weights 



Figure 2.18: Output Weightings for Design Narrow Band Disturbance 
Rejection 

changing the weighting on the controlled input u, but the weighting on y\ had 
to be so relaxed that the solution was to use a very low gain controller. Also, 
the weighting on y\ could be raised, but at the expense of having to decrease 
the controlled input weighting. The result shown here is a compromise of 
these two situations. 

The closed loop and open loop frequency response is shown for comparison 
in figure 2.19. 

A frequency response plot of the controller is shown in figure 2.20. 
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H-inf Design: Narrow Band Disturbance Rejection 
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Figure 2.19: Comparison of Open Loop and Closed Loop Transfer Functions 
for H 0 o Design Narrow Band Disturbance Rejection 
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Figure 2.20: Controller TYansfer Function for Design Narrow Band Dis- 
turbance Rejection 
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2.4 


Primary and Secondary Feedback vs Pri- 
mary Feedback Control 

This section explores the advantages of using primary and secondary feed- 
bank as opposed to only primary feedback. The primary feedback system is 
equivalent to setting K 2 equivalent to zero in figure 2.8. With the primary 
feedback system, the fed back quantity, yi, is the same quantity which the 
controller is trying to make zero. This presents a fundamental problem since 
the controller would “stop working” if yi were actually driven to zero. For in- 
stance, in the case of the disturbance d being a single sine wave at frequency 
10 Hz, the primary feedback system must have a greater gain in the feedback 
control in order to obtain the same level of performance as the primary plus 
secondary feedback system, because the two feedback system does not only 
rely on the measurement which the system is trying to drive to zero, y it but 
also a measurement correlated to the disturbance, namely y 2 . Following is 
an example of this. 

The following example uses the input and output weightings of the H 2 
design narrow band disturbance rejection example of the previous section. 
First, the plot of the closed loop frequency response from the disturbance d 
to the output j/i of both the primary and primary plus secondary feedback 
systems is shown in figure 2.21, to confirm that the performance levels are 
comparable. The control transfer functions are shown in figure 2.22. These 
controllers were designed for a narrow band disturbance acting at 10 Hz. 
The primary feedback controller is obviously at a much higher gain than the 
primary plus secondary feedback controller. The disadvantage of a high gain 
controller is that its performance is greatly effected by noise in the system. 


2.5 Conclusions for Part II 

A method for design of controllers for both wide and narrow band distur- 
bance rejection has been developed. This procedure allows the designer to 
incorporate a priori knowledge concerning the exogenous inputs, and to en- 
force design criteria in a systematic manner, thereby removing much of the 
“guess work” present in classical design methods. These designs utilized the 
H 2 and Hoc norms by using the Matlab Robust-Control Toolbox. 
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H2 Design: Narrow Band Disturbance Rejection: 1 Feedback vs 2 Feedback 



Figure 2.21: Comparison of Performance for the Primary (1) Feedback Sys- 
tem and the Primary Plus Secondary (2) Feedback System 
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Figure 2.22: Comparison of Controllers for the Primary (1) Feedback System 
and the Primary Plus Secondary (2) Feedback System 
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Also, the advantages of using a primary plus secondary feedback in the 
closed loop system for disturbance rejection was explored. The primary feed- 
back corresponds to the output which the control is attempting to control, 
while the secondary feedback corresponds to the output which is correlated 
to the disturbance. It was shown that a lower gain controller can be found 
when feeding back both primary and secondary outputs, as compared with 
feeding back only the primary output. 
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